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We review results from our monitoring observations of several lensed quasars performed in the optical, UV, and X-ray 
bands. Modeling of the multi-wavelength light curves provides constraints on the extent of the optical, UV, and X-ray 
emission regions. One of the important results of our analysis is that the optical sizes as inferred from the microlensing 
analysis are significantly larger than those predicted by the theoretical-thin-disk estimate. In a few cases we also constrain 
the slope of the size-wavelength relation. Our size constraints of the soft and hard X-ray emission regions of quasars 
indicate that in some objects of our sample the hard X-ray emission region is more compact than the soft and in others 
the soft emission region is smaller. This difference may be the result of the relative strengths of the disk-reflected (harder 
and extended) versus corona-direct (softer and compact) components in the quasars of our sample. Finally, we present 
the analysis of several strong microlensing events where we detect an evolution of the relativistic Fe line profile as the 
magnification caustic traverses the accretion disk. These caustic crossings are used to provide constraints on the innermost 
stable circular orbit (ISCO) radius and the accretion disk inclination angle of the black hole in quasar RX J1131 —1231. 
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1 Introduction 

Due to gravitational lensing the event horizon of a black 
hole casts a shadow that for a distant observer has a size of 
about lOfg (e.g., Falcke et al. 2000). The predicted angular 
sizes of these black-hole shadows for SgrA* and M87 based 
on their measured black-hole masses are ~ 50p,as and ^ 
35p,as, respectively. Direct imaging has been employed to 
image the accretion disk of SgrA* with submm Very Long 
Baseline Interferometry (Doeleman et al. 2008, 2011). 

Direct imaging of the environments of quasars using 
submm VLBI is not possible due to their large distances, 
however, we currently rely on indirect methods to infer the 
structure near black holes. Indirect mapping methods in¬ 
clude light travel time arguments, reverberation mapping 
of the broad line region (Blandford & McKee 1982; Peter¬ 
son 1993, Netzer & Peterson 1997), reverberation mapping 
of the Fe Ka emission region (Young & Reynolds 2000), 
reverberation between the X-ray and optical/UV continua 
(Shappee et al. 2014; Edelson et al. 2015), and microlensing 
of the continuum and line emission regions (e.g., Grieger et 
al. 1988 and 1991; Schneider, Ehlers & Ealco 1992; Gould 
& Gaudi 1997; Agol & Krolik 1999; Yonehara et al. 1999; 
Mineshige & Yonehara 1999; Chartas et al. 2002, 2009, 
2012; Popovic et al. 2003, 2006; Blackburne et al. 2006, 
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2011, 2014, 2015; Pooley et al. 2006, 2007; Kochanek et 
al. 2004, 2007; Jovanovic et al. 2008; Morgan et al. 2006, 
2008a, 2008b, 2010, 2012; Dai et al. 2010; Mosquera et al. 
2011, 2013; Chen et al. 2011, 2012; MacLeod et al. 2015). 

Macrolensing of a quasar into multiple images may oc¬ 
cur in cases where there is near alignment of the observer, an 
intervening galaxy and the background quasar. The surface 
mass density of the lensing galaxy needs to be above a criti¬ 
cal value before multiple images are produced (strong lens¬ 
ing). Typical separations of images produced in macrolens¬ 
ing range between 0.1 - 20 arcsec. 

Microlensing is the bending of light produced by the in¬ 
dividual stars in the lensing galaxy. Microlensing variability 
occurs when the complex pattern of caustics produced by 
stars in the lens moves across the source plane. 

The characteristic scale of these patterns is the Einstein 
radius of 


Re 


Dos 


- AG{M) Des 

DoeDos_ 


( 1 ) 


where (M) is the mean mass of the lensing stars, D is the 
angular diameter distance, and the subscripts L, S, and O 
refer to the lens, source, and observer, respectively. 

Microlensing will affect the images differently result¬ 
ing in uncorrelated variability between images. Emission 
regions with sizes significantly smaller than the projected 
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Einstein radius of the stars will be strongly affected by mi- 
crolensing whereas emission regions with sizes significantly 
larger will not be affected. The light curves of compact 
sources are therefore expected to show enhanced uncorre¬ 
lated variability compared to the light-curves of more ex¬ 
tended emission regions. 

Simulations of light-curves of images of gravitationally 
lensed quasars from caustic crossings are fit to observed 
light-curves and provide constraints on the size of the emis¬ 
sion regions, the microlens mass scale, and the mass frac¬ 
tion of the local surface density of the lens galaxy. The 
microlensing analysis includes the creation of many ran¬ 
dom realizations of the star fields near each image and the 
generation of magnification maps. This technique was first 
developed by Kochanek (2004; Kochanek et al. 2007) and 
has successfully been applied to several gravitational lens 
systems. Recent improvements of microlensing simulations 
that allow for movement of the stars between epochs also 
provide constraints on the inclination of the accretion disk 
and the direction of motion of the caustics (e.g., Poindexter 
& Kochanek 2010). 

The majority of the optical and UV continuum emis¬ 
sion of AGN is thought to originate from the accretion disk. 
The microlensing method therefore applied to the optical 
and UV light curves of images of lensed quasars places 
constraints on the sizes of the accretion disks at their re¬ 
spective rest-frame wavelengths. The majority of the X-ray 
continuum of quasars is dominated by non-thermal radia¬ 
tion of the X-ray corona with a smaller contribution from 
the disk. The accretion disk components is often referred 
to as the reflected component. Measurements using the mi¬ 
crolensing method to fit the X-ray light curves of the images 
of lensed quasars place constraints on the size of the X-ray 
hot corona. 

A promising technique for measuring the ISCO and spin 
of AGN relies on modeling the Fe Ka fluorescence lines 
originating from the inner parts of the disk (e.g., Fabian et 
al. 1989; Faor 1991; Reynolds & Nowak 2003). Modeling 
of the relativistic iron line can only be applied to a few rela¬ 
tively nearby Seyferts where the line is detectable at a high 
signal-to-noise level. The Fe Ka line is most Seyferts is 
typically very weak and constraining the spin and accretion 
disk parameters of Seyferts requires considerable observing 
time on XMM-Newton and Chandra. 

Here we provide for the first time a new technique based 
on microlensing that provides a robust upper limit on the 
size of the inner radius of the X-ray emitting region of an 
accretion disk. Whereas relativistically broadened Fe Ka 
lines detected in the spectra of unlensed AGN are produced 
from emission originating from a large range of azimuthal 
angles and radii, the microlensed Fe Ka lines in the spec¬ 
tra of lensed AGN are produced from a relatively smaller 
region on the disk that is magnified as a microlensing caus¬ 
tic crosses the disk. Simulations of microlensing caustics 
crossing the accretion disk show significant redshifts and 
blueshifts of the fluorescent Fe Ka line (e.g., Popovic et al. 


Table 1 Gravitationally Tensed Quasars Monitored 


Object 

Zs 

Zi 

Mbh 
(10« Mq) 

Line 

Q J0158-4325 

1.294 

0.317 

1.6 

(Mgll, a) 

HE 0435-1223 

1.689 

0.454 

5.0 

(CIV, a) 

SDSS0924-H0219 

1.524 

0.390 

2.8 

(Mgll, b) 

SDSS1004-H4112 

1.734 

0.680 

0.4 

(Mgll, b) 

QSO 1104-1805 

2.319 

0.729 

5.9 

(H^, c) 

PG 1115-1-080 

1.720 

0.310 

4.6 

(H/i, c) 

RX J1131-1231 

0.658 

0.295 

0.6 

(H/3, a) 

Q 2237-1-030 

1.690 

0.040 

12.0 

(H/l, c) 


Notes: 

a-Peng et al. 2006, b-Morgan et al. 2006, c-Assef et al. 2011 


2006). The evolution of the energy and shape of the Fe Ka 
line during a caustic crossing depends on the ISCO, spin, 
inclination angle of the disk and caustic angle. The extreme 
shifts are produced when the microlensing caustic is near 
the ISCO of the black hole. Measurements of the distribu¬ 
tion of the energy shifts of the Fe Ka line due to microlens¬ 
ing therefore provides a powerful and robust method for es¬ 
timating the ISCO, spin and inclination angle of the disk. 
In §2 we present a review of our multi-wavelength moni¬ 
toring campaign of gravitationally lensed quasars and sum¬ 
marise our major results. In §3 we present an estimate of 
the ISCO of RX J1I3I —1231 based on the measured dis¬ 
tribution of energy shifts of the Fe Ka line caused by mi¬ 
crolensing. Finally in §4 we present a summary of results 
obtained from our monitoring campaign of lensed quasars. 
Throughout this paper we adopt a flat A cosmology with 
Ho = 67 km s-^ Mpc-\ Oa = 0.69, and Om = 0.31, 
based on the Planck 2015 results (Planck Collaboration et 
al. 2015). 

2 Multiwavelength Monitoring of Lensed 
Quasars 

We are performing multiwavelength monitoring of several 
quasars listed in Table [T]with the main scientific goal of 
measuring the emission structure near the black holes in the 
optical, UV, and X-ray bands in order to test accretion disk 
models. The X-ray monitoring observations were performed 
with the Chandra X-ray Observatory. The Optical (B, R and 
I band) observations were made with the SMARTS Consor¬ 
tium 1.3m telescope in Chile. The UV observations were 
performed with the Hubble Space Telescope. Here we high¬ 
light several of the most recent results published by our mi¬ 
crolensing team. 

In Blackburne et al. 2015 we analyze the lightcurves 
of z = 2.32 quasar HE 1104—1805 using dynamical mi¬ 
crolensing magnification patterns. HE 1104—1805 has been 
observed at a variety of wavelengths ranging from the mid- 
infrared to X-ray for nearly 20 years. Dynamical microlens¬ 
ing analysis constrains the half-light radius of the accre¬ 
tion disk to log(ri/ 2 /cm) = 16.0lg gg , the half-light ra- 
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Fig.l X-ray half-light radii of quasars as determined 
from our microlensing analysis versus their black hole 
masses. 

dius of the X-ray emission region to have an upper limit 
of log(ri/ 2 /cm) = 15.33 (95% confidence), a low inclination 
angle is preferred statistically, the mean mass of the stars in 
the lensing galaxy (<M>) ranges between 0.1 and 0.4 M© 
and the slope of the size-wavelength relation oc A^, is 
^ = l.Od^o lg. The majority of the observed continuum X- 
ray emission is found to originate within ^ 30rg, assum¬ 
ing a black hole estimate of Mbh = 5.9 x 10® M 0 based 
on the width of the H/3 line (Assef et al. 2011). Based on 
this black hole mass estimate the gravitational radius of 
HE 1104-1805 is rg = 8.7 x lO^® cm. 

In MacLeod et al. 2015 we analyze the light-curves of 
2 : = 1.524 quasar SDSS 0924-1-0219 using static microlens¬ 
ing magnification patterns. SDSS J0924-I-0219 has been ob¬ 
served at a variety of wavelengths ranging from the near- 
infrared to X-ray. Our microlensing analysis in this system 
constrains the soft-X-ray, UV, and optical half-light radii to 
be 2 . 5 ^ 2 ” X10^^ cm, 8lf x lO^"^ cm, and - 5t^;5 x 10^® 
cm, respectively. Assuming the Mgll based black-hole es¬ 
timate of Mbh = 2.8 x 10® Mq the majority of the soft 
X-ray emission of SDSS 0924-1-0219 originates within ~ 
30 Tg. The gravitational radius of SDSS 09244-0219 is rg = 
4.12 X 10^® cm. 

In Dai et al. 2010 and Chartas et al. 2009 we analyze the 
light-curves of 2 = 0.658 quasar RX J1131 —1231. We find 
the X-ray and optical half-light radii to be 2.3 x 10®^ cm 
and 1.3 X 10®® cm, respectively. These sizes correspond to 
~ 26rg and ^ 147rg, respectively. 

An important result found in all microlensing studies is 
that optical sizes of quasar accretion disks as inferred from 
the microlensing analysis are significantly larger than those 
predicted by thin-disk theory. Specifically, measurements of 
the radius of the accretion disk at 2500A rest-frame indicate 
that the sizes obtained from microlensing measurements are 
2-3 times larger than the values predicted by thin-disk the¬ 
ory (e.g., Morgan et al. 2010, Mosquera et al. 2013). 

In Figure 1 we present the X-ray half-light radii of 
quasars from recent microlensing studies of lensed systems 
observed as part of our monitoring program (MacLeod et al. 
2015, Blackburne et al. 2014, 2015, Mosquera et al. 2013, 
Morgan et al. 2008, 2012, Dai et al. 2010, and Chartas et 
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Fig. 2 Evolution of the Fe Ka line possibly caused by the 
motion of a magnification caustic as it moves away from the 
center of the black hole. 

al. 2009). Included in Figure 1 are the uncertainties of the 
black-hole mass estimates and uncertainties in the size es¬ 
timates. The X-ray sizes of the quasars in our sample are 
found to be close to the sizes of their innermost stable circu¬ 
lar orbits. Assuming that most of the X-ray emission in the 
band detected originates from the hot X-ray corona, these 
results indicate that the corona is very compact and not ex¬ 
tended over a large portion of the accretion disk. 

3 Estimating the Inner Most Stable Circular 
Orbit Using Microlensing 

RX J1131 —1231 has been monitored 38 times over a pe¬ 
riod of 10 years with the Chandra X-ray Observatory. As 
reported in Chartas el al. 2012, redshifted and blueshifted 
Fe Ka lines have been detected in the spectra of the lensed 
images. 

In Figure 2 we show the evolution of the red and blue 
components of the Fe Ka line possibly caused by the mo¬ 
tion of a magnification caustic as it moves away from the 
center of the black hole. We interpret the shift of the Fe Ka 
line as resulting from general relativistic and special rela¬ 
tivistic Doppler effects. As shown in Figure 3, the two red- 
shifted iron lines in the Jan 1, 2007 observation are each 
detected at the >99% confidence and the iron lines in the 
Feb 13, 2007 observation are each marginally detected at 
the > 90% confidence level (Figure 3). 
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Fig. 3 The contours represent 68%, 90%, and 99% 
confidence intervals of the flux normalizations of the de¬ 
tected Fe Ka line in the spectra (Figure 2) of image B for 
the Jan 1, and Feb 13 observations of RX J1131 —1231. 



Fig. 4 Extremal shifts of the observed Fe Ka line energy 
for spin values ranging between 0.1 and 0.998 in increments 
of 0.1. Horizontal dashed lines represent the observed val¬ 
ues of the energy shifts g = EohJE^est of the most redshifted 
Fe Ka line components of the two epochs presented in Fig¬ 
ure 2. The observed shifts of the energy of the Fe Ka line 
can be interpreted as the motion of a magnification caustic 
over a distance of 16 ig between the two epochs. 
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Fig. 5 Distribution of the energy-shifts of the Fe Ka line 
in all images from all 38 epochs of RX J1131 —1231. Only 
cases where the iron line is detected at > 99% (panel a) and 
at > 90% (panel b) are shown. The vertical lines represent 
the extreme cut-offs of the distribution. These cut-offs pro¬ 
vide upper limits of the ISCO and inclination angle of the 
black hole. 


The extreme values of the normalized energy shifts g = 
£/obs/£'rest of the two lines can be used to roughly constrain 
the radius of the emitting material as shown in Figure 4. The 
generalized Doppler shifts g were calculated for the case of 
a spinning black hole (Kerr 1963) using the formalism de¬ 
scribed in Muller & Camenzind 2004 and Karas & Sochora 
2010 . 

For these calculations we assumed the radial and 
toroidal velocity components and v^^^ of the radiating 
plasma of the accretion disk in the Zero Angular Momen¬ 
tum Observer frame to be relatively small. For a given incli¬ 
nation angle i, the angle between the direction of the orbital 
velocity of the plasma and our line of sight is assumed 
to range between a minimum and maximum value that will 
depend on the inclination angle i, and the angle between the 
caustic direction of motion and the projection of our line- 
of-sight onto the disk plane. We assume that the most red- 
shifted line component, as shown in Figure 2, is produced 
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Fig. 6 Extremal shifts of the observed Fe Ka line energy 
for spin values ranging between 0.1 and 0.998 in increments 
of 0.1. Horizontal solid lines represent the observed values 
of the normalized energy shifts g = £’obs/£'rest of the most 
redshifted and blueshifted Fe Ka lines from all 38 epochs 
and all images of RX J1131 —1231. The extreme g values 
are for Fe Ka lines detected at > 99%. The inner radius of 
the accretion disk is constrained to be r < 8.5rg. 

by Fe Ka emission originating close to the ISCO. For r > 
r^s a Keplerian profile is assumed (equation 15 of Muller 
& Camenzind 2004), whereas for r < r^s constant specific 
angular momentum is assumed. 

The two redshifted iron lines detected in the Jan 1, 
2007 observation correspond to generalized Doppler shifts 
of Pmin ~ 0.7 and pmax 0.82 and the two lines detected 
in the Feb 13, 2007 observation correspond to pmin ^ 0.9 
and Pmax ~ 0.96. The detected energies of the iron lines 
and their evolution is consistent with a microlensing caustic 
moving by a distance of about 15rg between the two epochs 
that are separated by 44 days (observed frame). 

We can independently check this estimate by estimating 
the distance a caustic in this system is expected to travel 
in 44 days. Due to the combined motions of the observer, 
lens, source, and stars in the lens galaxy, the source moves 
relative to the magnification patterns at an effective velocity 
of Veil ^ 700 km s“^ (see Mosquera et al. 2011). Assuming 
a black-hole mass of ^ 6 x 10^ Mq (i.e., Sluse et al. 2012) 
the gravitational radius of RX J1131 —1231 is rg = 3.6 x 
10^^ cm. The distance Adc travelled by a caustic in Afobs = 
44 days is Adc = Veff Afobs/(l+'Z) ^ l^rg. This is consistent 
with the estimated value based on the shift of the Fe Ka 
lines between the two epochs. 

We have performed a systematic spectral analysis of all 
epochs and all images searching for the presence of the 
Fe Ka lines in the spectra. To determine the significance 
of the iron lines detected we produced confidence con¬ 
tours of the flux normalization of the lines as a function of 
their energy (see Figure 3). The spectral models used to fit 
the X-ray spectra of RX J1131 —1231 consist of a power law 


with Galactic and neutral intrinsic absorption and Gaussian 
lines to model the iron lines. For producing the confi¬ 
dence contours the only model parameter that was frozen in 
the spectral fit was the Galactic column density of Ah = 3.6 
X 10^'^ cm“^ (Dickey & Fockman 1990). 

We find that the iron line is detected in 58 out of the 152 
spectra (38 epochs x 4 images) at > 90% confidence and in 
18/152 spectra at > 99% confidence. There is considerable 
shift of the energy of the iron line due to microlensing. 

In Figure 5 we show the distribution of the energy-shifts 
of the Fe Ka line in all images. One important feature of 
this iron line energy-shift distribution in the significant cut¬ 
off of the distribution at rest-frame energies of ^ 3.9 keV 
and ^ 7.4 keV for iron lines detected at > 99% confidence. 
These cut-offs represent the most extremely redshifted and 
blueshifted Fe Ka lines. If we interpret the largest energy- 
shifts to be produced by X-ray emission originating close to 
the inner most stable circular orbit of the black hole we can 
provide upper limits on the size of the ISCO and inclination 
angle of RX J1131-1231. 

The extreme blueshift of the g distribution is sensitive to 
the inclination angle of the disk. Specifically, the measured 
generalized Doppler factor pmax = 115 (99% confidence) 
constrains the inclination angle to be > 55°. 

The extreme redshift of the g distribution is sensitive 
to the ISCO radius. Specifically, the measured general¬ 
ized Doppler factor pmin = 0.61 (99% confidence) con¬ 
strains 3.5rg < nsco < 9rg. For the measured generalized 
Doppler factor pmin = 0.53 (90% confidence) we find 3rg < 
DSCO < 7rg. 

In Figure 6 we plot the g values of the iron line for an 
accretion disk with an inclination angle of / = 60° and a 
caustic crossing angle of 9c = 20°. The observed extreme 
values of pmin = 0.6 and pmax =1.15 from all epochs and 
all images are consistent with an emitting region of an nsco 
< 9rg and an inclination angle of ^ 60°. 

A more detailed analysis will involve simulating the 
shifts of the Fe Ka line due to microlensing caustic cross¬ 
ings for a range of inclination angles, caustic crossing an¬ 
gles and black hole spin parameters. A comparison between 
the observed generalized Doppler shift g distribution with 
those simulated will provide more accurate constraints on 
the inclination angle, and the ISCO. The spin of the black 
hole can be inferred from the known relation between ISCO 
radius and spin. 

Additional monitoring observations with Chandra will 
provide more representative and complete distributions of 
the generalized Doppler shift g values of the Fe Ka line 
in RX J1131 —1231 and other lensed quasars. The g- 
distribution method provides a new and independent tech¬ 
nique of constraining the inclination angle, ISCO radius, 
and black hole spin of quasars. 
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4 Conclusions 

Multiwavelegth monitoring of lensed quasars provides a 
powerful technique for constraining the sizes of accretion 
disks and hot X-ray coronae of distant quasars. Simulta¬ 
neous modeling of the optical, UV and X-ray light curves 
places strong constraints on their relative sizes. Recent im¬ 
provements to this method include dynamical microlensing 
which provides additional constraints on the inclination an¬ 
gle of the disk. Several important results of our monitoring 
campaign of lensed quasars are: 

1. Optical sizes of accretion disks at 2500A inferred 
from our microlensing analysis are a factor of 2-3 larger 
than those predicted by thin disk theory. This result is 
consistent with recent measurements of continuum lags 
in nearby Seyferts NGC2617 (Shappee et al. 2014) and 
NGC5548 (Edelson et al. 2015). 

2. The size of the UV emission region is found to be a 
factor of about 10 on average smaller than the optical but 
larger than the X-ray emission region. 

2. The X-ray corona is found to be compact with a half¬ 
flight radius < 30rg (see Figure 1). 

4. The scaling between optical sizes at 2500A and black 
hole mass is consistent with thin disk theory. 

Future improvements include monitoring a larger vari¬ 
ety of quasars (for example, BAF quasars and radio loud 
quasars.) Several models predict that radio-loud quasars 
contain truncated disks and measurements of the sizes of 
radio-loud objects may address this issue. 

We introduced a promising new technique (g- 
distribution method) for measuring the inclination angle, 
ISCO, and spin of a black hole. The p-distribution method 
involves measurements of the distribution of the energy 
shifts of the relativistic iron line emitted from the accretion 
disk and a comparison of the measured ^-distribution with 
microlensing caustic simulations. The method has been 
applied to RX J1131 —1231 and initial results indicate 
that risco ^ 9rg and i > 60°. Further monitoring of 
RX J1131 —1231 and other lensed quasars will provide 
tighter constraints on the inclination angle, ISCO radius 
and spin of the black hole of distant quasars. 

Acknowledgements. We acknowledge financial support from 
NASA via the Smithsonian Institution grant G03-14110B. 
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